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INTRODUCTION 


I . 


Previous  work  done  in  connection  with  the  present  research  program 
has  utilized  two  modern  computational  methods  in  electromagnetics.  The 
first  technique  is  the  Geometrical  Theory  of  Diffraction  or  GTD,  which  is 
most  often  used  when  an  antenna  radiates  from  a structure  that  is  large 
in  terms  of  the  wavelength.  The  second  technique  is  the  Method  of  Moments 
which,  for  practical  reasons  based  on  computer  storage,  applies  best  to 
structures  that  are  not  large  in  terms  of  the  wavelength. 

A major  innovation  of  the  research  performed  in  connection  with  the 
prediction  of_satellite  antenna  radiation  patterns  lies  in  the  formal 
combination  of  the  GTD  and  the  moment  method  into  a single  computational 
method  known  as  the  Hybrid  Technique  [1].  This  technique  was  used  by 
Moore  and  Thiele  [2]  to  investigate  a turnstyle  type  antenna  on  a 24"  x 24" 
satel 1 i te  panel . 

The  computer  program  developed  in  connection  with  the  turnstyle 
antenna  investigation  was  limited  to  monopole  type  elements  and  pattern 
calculations  in  one  plane  only.  Subsequent  work  described  in  this  report 
has  oeen  concerned  with  the  general i zation  of  the  theory  and  the 
incorporation  of  that  theory  into  the  computer  program  so  that  it  may 
be  used  in  the  analysis  and  design  of  more  general  satellite  antenna 
systems  than  previously  possible. 

II.  THE  HYBRID  TECHNIQUE 

The  Hybrid  Technique  is  actually  an  extension  of  the  moment  method 
accomplished  by  modifying  the  impedance  matrix  to  include  the  GTD. 

A . Momen  t f^t h^^ 

Consider  a general  perfectly  conducting  body  on  which  there  is  a 
current  density  J.  In  order  to  satisfy  the  boundary  condition  everywhere 
on  the  body  surface,  the  total  tangential  electric  field  has  to  be  zero. 
That  is. 


= 0 


(2-1) 


where  eJ^^^  is  the  tangential  component  of  the  incident  electric  field 
due  to  a source  located  anywhere  outside  the  body  and  E^^^  is  the  tangential 
component  of  the  scattered  electric  field  radiated  by  the  current  J on  the 
conducting  body.  For  simplicity,  the  subscript  "tan"  is  dropped  and  it 
is  understood  that  only  the  tangential  electric  field  is  being  applied  to 
satisfy  the  boundary  condition  on  the  body  surface. 

To  apply  the  moment  method  to  find  the  unknown  response  J with  E^  as 
the  known  excitatign,  write  the  functional  operator  equation  [3,4]  to 
relate  the  J and  P as 


1 


{2-2) 


I (J)  t' 
op  ' 

whorp  Lop  r(?prosents  the  integral  functional  operator.  For  a unique 
solution,  L^p  must  have  a corresponding  inverse  operator,  l-op'^’ 

L^p'^  iC)  --  J . {2-2) 

In  addition,  an  inner  product,  obtained  by  integrating  J-E  over  the 
surface  under  consideration  to  obtain  the  reaction  J,F  [5],  is  defined 
to  satisfy  the  following  symmetric  inner  product  axioms: 

tj.|  + E'  ” ''BJ2»E  ' (2-^) 


J,E-  ^ E,J 


(2-5) 


J*,J  = 0 if  and  only  if  J =0 

■ n if  J 0 (2  6) 

where  » and  . are  scalars  and  * denotes  complex  conjugation. 

Next  we  expand  the  current  J on  the  surface  of  the  conducting  body 
in  a series  of  basis  functions  U,,  J^,  •••  J * defined  in  the  domain  of 
l-op-  ThPt  is, 

J = i 1„  J„  (2-n 


where  are  the  complex  coefficients  to  be  determined. 


Substituting  Equation  (2-7)  into  Equation  (2-2),  using  the  linearity 
of  the  operator,  and  forming  the  inner  product  with  a set  of  weight 
functions  IW,|,  •••,  W I in  the  domain  of  L , we  then  obtain 


N 

/. 

n=l 


1 -W  , 
n m 


Lop(^n) 


^m’ 


(2-8) 

I 


2 


Actually,  Equation  (2-8)  represents  the  nr"  row  of  equations  in  a system 
of  N such  equations  while  the  quantity  Lop(Jn)  represents  the  electric 
field  from  the  nth  basis  function  of  unit  amplitude.  In  the  usual  moment 
method  matrix  notation.  Equation  (2-8)  can  be  written  as 

[Z]  [I]  = [V]  (2-9) 


where  [Z],  [I]  and  [V]  are  called  the  generalized  impedance  matrix,  the 
generalized  current  matrix  and  generalized  voltage  matrix,  respectively. 
The  elements  of  [Z]  are  denoted  by 


Z 

mn 


(2-10) 


^ In 

To  include  the  diffracted  field  contribution  from  the  n'^  basis 
function  J^,  a new  matrix  element  can  be  written  as 


V 

mn 


L 

op 


(J  ) + c L (J  )> 
' n'  op'  n' 


or,  by  axioms  (2-4)  and  (2-5), 


(2-11) 


'mn 


= <W 


+ 


(2-12) 


where  c Lop(Jn)  f'epresents  the  contribution  from  GTD  and  c is  a complex 
quantity.  Equation  (2-12)  may  be  further  simplified  as 


V = Z + Z^ 
mn  mn  mn 


(2-13) 


where  Z,^p  = -Wpi,  c Lop(Jn)^-  With  the  new  matrix  element.  Equation  (2-9) 
can  be  rewritten  as 


[Z-]  [r]  = [V] 


(2-14) 


where  [I']  is  the  current  distribution  on  the  conducting  body  surface, 
including  the  contribution  from  GTD. 
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B.  [)  1 f f r'fic tion  Theory 

The  problem  of  straight  edge  diffraction  by  a perfectly-conducting 
wedge  was  first  solved  by  Sommerfeld  [6],  Later,  Pauli  [7]  introduced 
the  Vp,  function  as  a practical  formulation  in  the  solution  for  the  fields 
diffracted  by  a perfectly-conducting  wedge  of  finite  angle.  More 
recently,  Hutchins  and  Kouyoumjian  presented  a formulation  which  yields 
better  accuracy  than  Pauli's  solution  in  the  transition  regions  near  the 
incident  and  reflected  shadow  boundaries,  particularly  when  r^ x . Their 
formulation  can  be  written  as 

V,^(L,i  ,n)  = I_^(L,I',n)  + l^^(L,f,n)  (2-lb) 

where 


I 


(L,r  ,n) 


p-j(kL+i./4) 

jr)J^ 


^ ('Tn  ) 


eJl^La 


and 


[higher  order  terms] 


a = 1 + cos(i-2nTiN) 


(2-16) 


(2-17) 


The  higher  order  terms  are  negligible  for  large  kL  and  n is  defined  from 
the  wedge  angle  WA  = (2-n)n  as  shown  in  Figure  2-1.  Moreover,  N is  a 
positive  or  negative  integer  or  zero  whichever  most  nearly  satisfies  the 
equation 


2nnN-p,  = -n  for  I (2-18) 

- Tt 

and 

2nu  N-p  = -n  for  I,  . (2-19) 

+ 7T 


The  distance  parameter  L and  H are  defined  subsequently. 


4 


<■9- 


Considi'r  fiqur('  ?-l  with  source  fieiil  C (s)  fron  a source  located 
at  point  s'(,',:',z').  The  source  can  be  either  el r>ctric  or  magnetic 
giving  rise  to  either  plane,  cylindrical,  conical  or  spherical  wave 
incidence  on  the  wedge  tip.  The  diffracted  vector  field  at  the  point 
s(.  ,:,7)  can  he  written  in  terms  of  a dyadic  diffraction  coefficient. 
Kouyounij  lan  and  Pathak  [S]  have  shown  that  thf’  diffracted  field  can  be 
written  compactly  if  they  are  expressed  in  terms  of  a ray  fixed 
coordinate  system  centered  at  the  point  or  points  of  diffraction  Q[;. 

The  rel ationshi ps  between  the  orthogonal  unit  vectors  associated  with 
these  coordinates  (s' ;s,f^,; ) are  given  by 

I = - s'  (2-20) 


1 = . ■ X (2-21) 


s = .'o  X f (2-22) 

where  I is  the  incident  direction  unit  vector,  and  s is  the  diffraction 
direction  unit  vector.  In  this  coordinate  system,  the  diffracted  field 
may  he  wri tten  as 


E^(s) 

e'(Q^) 

• 5(s,l)  A(s) 

(2-23) 

or,  in  terms 

of  the 

Vg  function  as 

o'  e;(Qj) 

(2  24) 

|o 

-vj 

0 

where 

■ 

e"^(s)  • 

(2-25) 

e‘'(s)  = 

E^(s)  • 

(2-26) 

1 


and 

Vg  = V^(L,i-',n)  ? Vg(L,f/,n) 


(?-^7) 


I 


I 

( 

I 

I 

1 

I 


Tt».'  minus  siqn  of  the  Vg  function  is  used  for  the  soft  case  (L  field 
parallel  to  the  edqe)  and  the  positive  siqn  is  for  the  hard  case 
((  field  perpendicular  to  the  edge).  Ttie  angular  relations  are  qiven 
by  ^ (pli'  where  the  tenii  is  associated  with  the  incident  field 

and  tl)e  t)  term  witti  the  reflected  field.  The  ray  divergence  factor  or 
spatial  attenuation  factor  A(s)  is  given  by 


A(s)  = 


J_ 

* s 


1 ^ V- 

i \ sTs+sT 


and  L is  given  by 
r 


. 2 

s sin  b. 


J 

. p ' +p 


s s sin 
s'+’s'' 


for  plane,  cylindrical  and 
conical  wave  incidence 


for  spherical  wave 
incidence 


for  plane  wave 
incidence 

for  cylindrical  wave 
incidence 

for  conical  and  spherical 
wave  incidence 


(2-28) 


2-29) 


It  is  the  introduction  of  the  distance  parameter  L by  Pathak  and 
Kouyoumjian  that  permits  near-zone  to  near-zone  diffraction  to  be  obtained 
accurately  for  kL>l . 

III.  APPLICATION  or  THL  HYBRID  TECHNIQUE 
Examples 

A.  Monopole  near  a Conducting  Wedge 

Consider  the  conical  problem  of  a iiionopole  near  a perfectly- 
conducting  wedge  as  shown  in  Figure  3-la.  Assume  that  the  monopole  is 
a thin  wire  with  radius,  R‘-^A  and  its  axis  coincides  with  the  coordinate 
z-axis.  To  find  the  current  distribution  on  the  monopole,  first  apply 
the  boundary  condition  on  the  wire  surface.  That  is, 

- r^(O.z)  = r (0,z)  (3-1) 

or 

L^p(J)  = - E^(0,z)  = eJ(0,z)  (3-2) 

where  E^(0,z)  may  be  expressed  as  [9] 
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Figiit'e  i-1.  (d)  Ilonopolfi  near  a porfectly  conducting  wedge 

(L))  The  direct  und  diffracted  fields  caused 
by  ttie  n^l'  source  segment  incident  upon 
the  mt^'  oliservation  segment. 


I 


l'^(U.z)  j J(z')  IP  [fl-jkr)(2r^-3K^)  + (kkr)^)]dz' 

' (3-3) 

where  z'  dnd  z denote  ttie  source  point  and  observation  point,  respectively, 
and 


r = + iz~z' 


(3-4) 


Therefore.  can  be  defined  as 


-Of) 


-A  r — 

Me 


G(r,r' ) dz' 


(3-b) 


where 


G(r,r')  = 


■ j k r 


[(l+jkr)(2r^-3R^)+(kRr)^] 


(3-b) 


Then  the  current  distribution  J(z')  may  be  e^anded_by  a finj_te  set  of 
piecewise  sinusoidal  basis  functions  [10],  { F-]  ( z ' ) , f2(  z ' ) • • ■ Ffj(  z ' ) i as 
shown  in  F inure  3-2.  That  is. 


fj  _ 

j(z')  = y I F (z‘) 

' n n 

n -1 


(3-7) 


■Jtiere  is  a complex  quantity  and 

sinh  k (z’-z^_.|)  sinh  k (z^^^^-z’)  , 

* f/ ^ ^ si nh  k (z  -z  ,T  ^ ^ sinh  k (z  ',-z  ) ^ 
n n-1  n+1  n 

i is  the  complex  propagation  constant  of  homogeneous  medium; 


(3-8) 


k = jw'ML 


(3-9) 


Tiius , equation  \3~2)  can  be  written  as 

: I i-  (r  (z‘ ))  = - r^(o,z)  . 

M op  n 
rr  1 


(3-10) 
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figure  3-2.  Piecewise  sinusoidal  basis  functions. 


Next,  witti  a set  of  weighing  functions  defined  as 


IT  = F (z' ) 

m ir 


(3-11) 


the  riatrix  element  in  Equation  (2-13)  can  be  calculated  by 


Z 

mn 


<F 

m 


(z’), 


L 

op 


(F  (z')> 
n ' 


(3-12) 


which  is  physically  interpreted  as  the  direct  field  contribution  at 
region  tn  due  to  a unit  of  piecewise  sinusoidal  current  Fn(z')  defined  at 
region  n as  shown  in  Figure  3-lb.  Also,  the  matrix  element  in 
Equation  (2-13)  that  accounts  for  the  diffracted  field  from  the  edge  of 
the  wedge  due  to  the  same  basis  function  Fp(z’)  can  be  written  as 


<F  (z 
m' 


). 


(3-13) 


where  ^(s),  the  diffracted  field  arrivino  at  renion  m,  is  obtained 
using  GTD  and  the  integral  is  over  the  length  of  region  m.  That  is, 

(3-H) 

0 

where  E^((J[^)  is  the  incident  field_component  perpendicular  to  the  edge 
due  to  the  current  basis  function  F^(z')  and  « is  defined  as  shown  in 
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I 

I 


Fiqure  3-lb.  In  this  particular  caso,  Ii|^|(Qt),  the  incident  field  con- 
ponent  parallel  to  the  edge  equals  zero,  and  the  incident  angle  is 
- -^ri.  Moreover,  for  sperical  wave  incident  A(s)  and  1.  are  defined 
by  Lquation  (2-28)  and  tquation  (2-29),  respectively.  Therefore, 
tquation  (3-14)  can  further  be  simplified  as 


I^(s) 


-V, 


jk[ 

e 


s ‘ s 
s ' +s 


(3-15) 


iilith  Z,!,„  ^ ^nr/^mn  t-^lculated,  the  unl.nown  current  distribution  on  the 
i.ionopole  can  then  be  obtained  by  solvinn  the  matrix  equation 

[!■]  = U'T^  [VJ  . (3-16) 


I 

I 

I 

I 
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B . Moriopoles  on  a Rectangular  Ground  Plane 

In  general,  given  the  locations  of  a source,  an  observation  point 
and  the  edge  of  a wedge,  it  is  necessary  to  determine  if  the  given  edge 
has  a diffraction  point  that  will  give  a diffraction  contribution  in 
the  observation  direction.  The  equations  derived  below  can  be  referred 
to  [11,12]. 

from  Figure  3-3,  the  position  vector  M of  a diffraction  point  can 
be  obtained  as  follows: 


RQ  = C?  - 

.[(^  - ^)-EJiE 

(3-17) 

U = (RO  - 

i^)-E 

(3-18) 

II 

UI 

(3-19) 

G = |R0  - 

RP 

(3-20) 

H = RT  - 

(3-21) 

RD  = RP  - 

(lili_)  L 

(3-22) 

Diffraction  occurs  only  if  the  following  inequality  is  satisfied. 


. L < 5-L  < . L 

,rr  - - rsi 


(3-23) 
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rigure  3-3.  Geometry  for  the  calculation  of  the  position 
vector  RD  of  a diffraction  point  Q^. 

After  RD  is  calculated,  the  same  procedure  as  in  Example  A will  give  the 
current  distributions  on  monopoles  on  the  ground  plane  in  the  presence 
of  the  diffracting  edges.  To  determine  the  far  field  pattern  including 
the  four  edges  of  diffraction,  it  is  necessary  to  define  the  incident  and 
reflection  shadow  boundaries  as  shown  in  Figure  3-4. 

From  Figure  3-5a,  it  is  obvious  that  only  if  the  sign  of 
n ■ - TT)  and  n-D  (3-24) 


are  the  same  and  the  intersection  point  is  not  on  the  rectangular  plate, 
then  tiie  observation  point  is  in  the  direct  incidence  region.  If  ^ 
represents  the  location  of  the  image  of  a source  and  if  the  intersection 
point  is  on  the  rectangular  plate,  then  the  opposite  sign  of  Equation  (3-24) 
indicates  that  the  observation  point  is  within  the  reflection  region  as 
shown  in  Figure  3-5b.  The  location  of  intersection  point  is  defined  as 
follows  (see  Figure  3-5c). 


= - (^  - rn-n 

D-n 


(3-25) 
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Figure  3-b(d).  Geometry  for  finding  direct  incidence  region. 


A 


IMAGE 


(b) 

Figure  3-b(b).  Geometry  for  finding  reflection  region. 
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SOURCE 


(C) 


Figure  3-6(c).  Geonietry  for  finding  the  location 
of  intersection  point. 


RT  = RS  + D . 

To  deteniiine  if  RT,  the  interception  point,  is  on  the  ground  plane 
the  following  equation  can  be  tested. 


n 


< Ti  not  on  the  plate 
> T on  the  plate 


where 


-1 


[((.„  - RT)  X - RT)] 


RTiT.n  1 


Figure 


Figure  3-6 


(o) 

3-G(b).  The  vector  relation  between  the  intersection 
point  and  one  ground  plane  corner. 

C3  C2 


(b) 

(b).  When  the  intersection  point  is  on  the  ground  plane, 

4 


> 71 


1 n=l 


figure  3-6(c).  When  the  intersection  point  is  not  on  the  ground  plane 


To  include  the  diffraction  by  the  four  edges  to  the  total  far  field 
pattern,  recjuiros  another  set  of  etjuations  to  deternine  the  diffraction 
points  on  tlie  four  edges,  sitice  in  far  field  case,  only  the  observation 
direction  is  given  and  it  is  assumed  that  r-<^. 


Figure  3-7.  Geometry  for  finding  the  position  vector  TO  of  a 
diffraction  point  when  the  far  field  observation 
direction  D is  given. 


From  Figure  3-7,  it  is  easy  to  show  that 


Q = 

(TO  - rr)  • E 

(3-29) 

M = 

1 (TO  - TO)  - Q E 1 

(3-30) 

cot 

D-E 

(3-31) 

^0  — - . - 2 
\ 1 - (D-E)'^ 

^ = rr  + g C 

TO  = ^ + (H  cot  L 


(3-32) 

(3-33) 


wherp  RD  is  the  location  of  diffraction  point  v;hich  nust  satisfy  the 
followinij  inequality  to  ensure  that  diffraction  occurs: 


-L  < D-L  < .[ 

Cl  - RS|  lC7  - RS  . 


(3-34) 


That  is,  the  diffraction  points  must  be  defined  on  the  edqes. 

The  total  far  field  pattern  is  calculated  by  sumninq  the  direct, 
reflected  and  diffracted  field.  The  direct  and  reflected  field  may  be 
obtained  via  the  moment  method  while  the  diffracted  field  nay  be 
determined  by  the  GTD.  As  mentioned  in  Section  II,  the  diffracted  field 
can  be  expressed  compactly  in  tenas  of  a fixed  ray  coordinate  system 
centered  at  the  point  Q^.  The  three  orthogonal  unit  vectors  of  the 
coordinate  system  as  shown  in  Figure  3-8  are  E the  unit  vector  along  the 


A 


I 

I 

I 

§ 

I 

I 

i 

1 

I 


i 


Figure  3-8.  Fixed  ray  coordinate  system  centered 
at  the  diffraction  point  Q£. 


I 

I 


1 


18 


oclqe,  n thp  non:ial  to  the  surface,  and.B  = n x [.  The  incident  unit 
vector  1,  the  diffraction  unit  vector  D,  dqi  r' 

represented  as  follov/s: 

- I = cos  p'  L + sin  c‘  cos  j'B  + smd'sin'j.'n  (-3-3B) 

0 0 0 

IJ  = cos  L + sin  cosi^  b + sin  b^  sin?  n (3-36) 


?•  = tan-1 

c-I-B/ 

(3-37) 

r \ 

. - 1 D • n 

? = tan  ' T— ^ 

■ J)'B/ 

(3-38) 

?'  = cos?'  n - sin?'  B 

(3-39) 

? = cos?  n - sin?  B 

(3-40) 

b'  = ?'  X I 

0 

(3-41) 

**  = .•  V n 

"o' 

(3-42) 

Thereupon,  the  diffracted  field  can  be  expressed  as 

7^  r ^ r ^ 

L = ? -►  E,,  B^ 

(3-43) 

where 


-V, 


; c;(o 


-J  L- 


rind 


= E’(Qp)  • 


I i , r 
i e;,(Qb) 


■'^bJ 


s’inB, 


A(s) 


(3-44) 


(3-45) 

(3-46) 
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wnero  l/'(Q[;)  is  the  incident  field  at  diffraction  point  Q[;  due  to  the 
current  distribution  on  the  monopoles.  For  spherical  wave  incident  and 

(3-44)  can  be  written  as 

jks'sin^c  jk^-D 
s'  e ° e (3-47) 


the  phase  is  referenced  to 

(x,y,z)  coordinate  systen  of  the  rectangular  plate  as  shown  in  Fiqure 
3-9,  the  phase  term  is  added. 


A 


Figure  3-9.  The  phase  reference  of  the  far  field  pattern  to 
the  coordinate  system  of  the  finite  ground  plane. 


The  diffracted  field  just  calculated  may  be  considered  as  arising 
from  a current  basis  function  InFn(*’)  defined  on  the  n^*^  region 
(Figure  3-9)  of  the  monopoles.  Therefore,  the  total  diffracted  field 


for  far  field  calculation  (s-*^),  Equation 


f^?,(s) 


L_L'i(s)_ 


! 

I 0 


-V 


B-J 


e1(qJ 


iU 


where  lim 


-jks 


is  suppressed  and,  since 


20 


ts  _the  sunma_ti on  of  the  contributions  fron  all  the  current  basis  functions 
•••  includinq  all  the  effects  from  the  four 

edqes  (excludinq  the  2nd  order  diffraction  and  the  corner  di f fraction) . 

In  general  , the  total  far  field  can  be  written  as 


L 


total 


N 

I 


n=l 


n 


r • r* 

n n 


(3-48) 


— 'j 

where  fp,  Lp  and  (,p  represent  respectively  the  direct  incident  field, 
refh!cted  field  and  diffracted  field  due  to  the  current  basis  function 
In  rn('")  defined  on  the  region  n of  the  monopoles. 

1 . Jingle  iionopole  on  a rectangular  plane 

r.onsider  a quarter  wavelength  monopole  excited  with  1 volt  generator 
on  a ground  plane  as  shown  in  Figure  3-10.  In  Figure  3-11,  the  calculated 
a-component  of  the  far  field  pattern  including  diffraction  by  the  4 edges 
is  coiiipared  with  that  obtained  by  Burnside's  GTD  computer  program  [13] 
for  the  following  cases: 


a) 

, f=900  MHz 

(Figure 

3-lla,b) 

b) 

Xi=><2=yi=y2=H" 

, f=800  MHz 

(Figure 

3-llc,d) 

c) 

, f=900  MHz 

(Figure 

3-lle,f) 

and 

d) 

, f=1000  MHz 

( F i gure 

3-llg,h) 

where  the  pattern  is  calculated  with  the  observation  direction  defined 
by  if  0”,  180"  and  0'’^U;^180".  In  Figure  3-llg,h,  the  discontinuity 
appearing  around  the  region,  0=90°  and  270°,  is  due  to  the  second  order 
diffraction  from  the  edges  which  is  not  included  in  the  calculations. 

The  reason  for  making  tiie  comparison  with  Burnside's  pure  GTD 
program  is  to  provide  an  independent  check  on  the  validity  of  our 
computer  program. 

Burnside's  program  uses  infinitesimally  short  current  sources  which 
have  amplitiides  and  locations  so  as  to  simulate  a quarter  wavelength 
monopole  witii  a sinusoidal  current  distribution. 
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OBSERVATION 

direction 


I VOLT  GENERATOR 


(b) 

Figure  3-10.  Geometry  for  one  inonopole  on  a finite  ground  plane. 
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0 io  «0  DO  so  too  1?0  1X0  ISO  ISO  ?00  720  7U3 


PfiTTERN  BNGLE.  THETfl  (DEGREES) 


MONOPflLE  TILTS  RT  90.0  DEGREE 
K(3NOPOLE  LENGTH-  D.35  WAVELENGTHS 
FREQUENCr-  900  MHZ 

O PATTERN  including  4 EDGES  OIFFf.RCTION 

Figure  3-1 1(a).  Ee  pattern  of  a stub  calculated  by  the  hybrid  technique. 


0 20  40  60  00  100  120  140  160  ISO 

PATTERN  ANGLE  6 

Figure  3-1 1(b).  Le  pattern  of  a stub  calculated  by  Burnside's  GTD  program 
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PATTERN  ANCLE.  THETA  (DECREES) 


I 


NONOPOLE  TILTS  AT  90.0  DEGREE 
HONOPOLE  LENGTH-  0.2S  HAVELENGTHS 
FREQUENCY-  600  MHZ 

C3PRTTERN  INCLUDING  Q EDGES  0IFFRRCTK5N 


Figure  3-1 1(c).  E,  pattern  of  a stub  calculated  by  the  hybrid  technique 


Figure  3-n(d).  pattern  of  a stub  calculated  by  Burnside's  GTD  program. 


HONOPOLE  TILTS  AT  90.0  DEGREE 
MONOPOLE  LENGTH-  0.25  WAVELENGTHS 
fREQUENCY-  900  MHI 

a PATTERN  INCLUDING  U EDGES  DIFFRACTION 


Figure  3-1 1(e).  pattern  of  a stub  calculated  by  the  hybrid  technique. 
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igure  3-ll(f).  E,  pattern  of  a stub  calculated  by  Burnside  s GTD  program. 


RELATIVE  SIGNAL  LEVEL  CCS) 


I 

I 
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MONOPOLE  TILTS  RT  90.0  DEGREE 
MONOPOLE  LENGTH-  0.25  WAVELENGTHS 
FREQUENCT-  1000  MHZ 

n PATTERN  INCLUDING  U EDGES  DIFFRACTION 


rifjure  3-ll(g).  pattern  of  a stub  calculated  by  the  hybrid  technique. 
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GTD  prog^an. 


1 1 . Two  tilted  nonopoles  on  a rectangular  plate 

Considor  two  quarter  wavelenqth  I'lonopoles  tilted  at  an  angle  of 
on  a ground  plane  excited  with  generators  as  shown  in  figure  3-12. 

The  i) -component  of  the  far  field  pattern  including  diffraction  by  the 
4 edges  is  compared  with  that  obtained  by  fiurnside's  GTD  computer  program 
for  the  following  cases: 

a)  x^-y^=8‘',  , f=1000  flHz  (Figure  3-13a,b) 

b)  x^=x^=y^=y2=l4"  , f=900  MHz  (figure  3-13c,d) 

c)  x^=ypl2",  X2=y2=16",  f=800  MHz  (Figure  3-13e,f) 

wiiere  the  observation  direction  is  given  by  1>=0°,  180°  and  0°^fi^l80°. 

It  is  noted  that  in  both  cases  i and  li,  the  i).- component  of  the 
far  field  pattern  observed  in  the  plane  l/=0°,  180"  and  0°^6M80°  is 
zero  since  the  antennas  in  these  configurations  have  no  contribution  to 
((.-component  in  that  particular  observation  direction.  In  order  to  see 
the  (^-component  of  the  far  field  pattern  of  two  tilted  antennas  as 
mentioned  in  case  ii,  the  observation  plane  (|>=90°,  270"  and  0°^6M80° 
is  used.  The  result  is  shown  in  Figure  3-14  for  the  cases: 

a)  x^=y^=8",  x^=y^=20"  . f=1000  MHz  (Figure  3-14a) 

b)  x^=y^=12",  X2=y2°T6">  f=800  MHz  (Figure  3-14b) 


C.  Helical  Antenna  on  a Ground  Plane 


The  helical  antenna  [14]  as  shown  in  Figure  3-15a  is 
the  following  symbols: 


described  by 


D = diameter  of  hel ix 

C = circumference  of  helix 

S = spacing  between  turns 

a = pitch  angle 

L = length  of  1 turn 

N = number  of  turns 

A = axial  length 

d = diameter  of  helix  conductor. 
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A 

z 


I 


I 

$•0* 

f i'jurt,"  J-IZ.  (d)  J-diiiioniional  (jt'niiietry  of  two  tilted  i;ionopolf?s  on  d 
finite  ground  pKine. 

([))  2-dimens  i oridl  gooiiietry  of  two  tilted  r.onopoles  on  d 
finite  ground  pi  due. 


1 1 1 1 T 1 1 1 1 1 1 1 1 1 1 p I r 

0 to  00  to  (0  too  120  lliO  too  1*0  200  220  2<t0  200  200  000  r'-'O  SOO  MO 

PflTTEBN  ANGLE.  THETA  (DEGREES) 


MONOPOLES  TILT  20.0  DEGREE 
M0N(5P0LE  LENGTH-  0.25  WAVELENGTHS 
PREQUENCY-  1000  MHZ 

GRP  BETWEEN  ANTENNAS  TIPS-  0.2  INCHES 
O PATTERN  INCLUDING  H EDGES  DIFFRACTION 


figure  3-13(a).  pattern  calculated  by  the  hybrid  technique. 


32 


Figure  3-13(b).  E pattern  calcul 


PftTTERN  ANGLE.  THETA  (DEGREES) 


MONOPOLES  TILT  AT  20.0  DEGREE 
HONOPOLE  LENGTH-  0.25  WAVELENGTHS 
FREQUENCY-  900  MHZ 

GAP  BETWEEN  ANTENNAS  TIPS-  0.2  INCHES 
ID  PATTERN  INCLUDING  U EDGES  DIFFRACTION 


Figure  3-13(c).  L pattern  calculated  by  the  hybrid  technique 


igure  3-13(,d). 


PATTERN  ANGLE.  THETA  tOEGREES) 
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HONOPOLES  Tl'.r  AT  20.0  DEGREE 
MaNOPOLE  LENGTH-  0.25  WAVELENGTHS 
FREQUENCY-  800  MHZ 

GAP  BETWEEN  ANTENNAS  TIPS-  0.2  INCHES 
D PATTERN  INCLUOING  <4  EDGES  DIFFRACTION 


Fiqure  3-13(e).  pattern  calculated  hy  the  hybrid  technique 


gure  3-13(f).  E pattern  calculated  by  Burnside's  GTD 


RELATIVE 
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T r r r 1 1 1 1 1 1 1 1 1 i t 

to  too  120  mo  too  leo  ?oo  zao  2uo  2so  2M  soo  S2o  no  mo 

PflTTCRN  ftNiSLe.  THETfl  (OEGREtSi 


MONflPflLES  U1.T  RT  20.0  DEGREE 
MONOPQLE  LENGTH"  0.25  HfiVELENGTHS 
FREQUENCY"  1000  MHZ 

GRP  BETWEEN  flNTENNRS  T1P3=  0.2  INCHES 
O PRTTERN  INCLUDING  >1  EDGES  DIFFRACTION 


I'igure  3-14(d) . 


L pattern  calculated  by  the  hyorid  technique. 

1' 
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PATTERN  RNbLE.  THETA  (DEGREES) 


-•e-  ISO* 


MONOPOLES  TIL-^AT  30.0  DEGREE 
MONOPOLE  LENGTH-  0.25  WAVELENGTHS 
FREQUENCY-  000  MHZ 
CAP  BETWEEN  ANTENNAS  TIPS-  0.2  INCHES 
O PATTERN  INCLUDING  U EDGES  DIFFRACTION 


Figure  3-14(b).  E,  pattern  calculated  by  the  hybrid  technique 


SURFACE  OF  IMAGINARY 
HELIX  CYLINDER  d 


rifjure  3-lb.  (a)  Helix  and  associated  dimension 

(b)  Relation  between  circumference, 
spacing,  turn  length,  and  pitch 
angle  of  one  turn  of  a helix. 
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If  1 turn  of  a circular  helix  is  stretched  out  straiqht  on  a flat  plane, 
the  relation  between  the  spacinn  5,  ci rcunference  C,  turn  length  L and 
pitch  angle  u are  as  illustrated  by  the  triangle  in  Finure  3-15b. 

For  the  work  presented  in  this  report,  a tiexagon  is  used  to 
approximate  the  circular  helical  antenna  which  has  the  teminal  ar- 
rangement as  shown  in  Figure  3-16a. 

In  general,  although  there  is  an  infinite  variety  of  radiating  far 
field  patterns  associated  with  a helical  antenna,  two  kinds  are  of 
particular  interest.  One  is  the  axial  or  beam  mode  of  radiation,  and  the 
other  is  the  normal  mode  of  radiation.  For  the  axial  mode  radiation, 
there  is  a field  maximum  in  the  direction  of  the  axis  of  the  helical 
antenna  and  for  tFie  nomial  mode  radiation,  the  field  maximum  occurs 
normal  to  the  helix  axis. 

To  obtain  the  far  field  pattern  of  a helical  antenna  on  a finite 
ground  plane  the  Hybrid  technique  is  applied.  The  far  field  patterns 
of  a 2-turn  fielix  radiating  in  an  axial  mode  on  an  infinite  ground  plane 
are  shown  in  Figure  3-17a  and  b;  whereas,  the  far  field  patterns  of  the 
same  helix  on  a finite  ground  plane  with  x-|=X2=y]=y2=0. 5>.  are  shown  in 
Fi(]ure  3-17c  and  d.  The  pitch  annle  and  circumference  of  the  radiating 
helix  are  13.965''’  and  1.0053  wavelencjth,  respectively.  The  lack  of 
complete  symmetry  is  due  to  the  physical  connection  between  the  helix 
itself  and  the  feed  point  which  renders  the  antenna  geometry  non- 
symmetrical . 

For  tne  normal  mode  radiation.  Figure  3-18a  and  Figure  3-18b  show 
the  far  field  patterns  of  a one-turn  helix  on  an  infinite  ground  plane  and 
a finite  ground  plane  (x-|=X2=yi=y2=0.5X) , respectively,  where  a=9.0341° 
and  c^O.ItiA.  The  discontinuity  appearing  around  b=90"  and  270°  is  due  to 
the  second  order  diffraction  from  the  edges  which  is  not  included  in  the 
calculations.  It  is  noted  that  only  the  F-e, -component  is  shown  in 
Figure  3-18.  Since  the  F^-component  is  practically  zero  compared  with 
F^-component  when  observed  in  a plane  with  4i=90°,  270°  and  0°^P^180°. 

Furthermore,  the  terminal  impedances  of  a 5-turn  helix  on  an 
infinite  ground  plane  and  on  a finite  ground  plane  with  xi=xiJ=yi=y2=0.5X 
are  calculated  witFi  different  frequencies  and  are  shown  in  Figure  3-19a 
and  b.  The  helix  is  approximated  by  hexagons  with  a=18°  and 
0.66X,'.c',  1 .2X . The  tenninal  impedance  is  highly  sensitive  to  changes 
in  frequency  for  c<2/3X,  but  is  nearly  as  constant  with  frequency  for 
3/4<C‘4/3  within  which  the  lielix  is  in  an  axial  mode.  The  spiral  form 
of  the  impedance  can  be  in  general  compared  to  the  impedance  measure- 
ment obtained  by  Kraus  [14]  as  shown  in  Figure  3-20.  The  disagreement 
between  Figure  3-19  and  Figure  3-20  is  due  to  the  sensitivity  of  the 
feed  point  modeling.  For  this  work,  a delta  gap  model  generator  is 
used  and  no  attempt  is  made  to  model  the  tapered  feed  in  the  vicinity 
of  the  feed  point  itself. 
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0 = 180° 


N = 2 0° 


Fiijure  3-17(a).  The  calculated  Ee  pattern  (not  in  dB)  of  a helical 
antenna  on  an  infinite  ground  plane  radiating  in 
the  axial  mode.  The  observation  plane  is  defined 
by  t=90°,  270°  and  0°<e<180°. 
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d =0.002X 

a = 13.965° 

C = I.0053X 
N = 2 

Figure  3-17(b).  The  calculated  pattern  (not  in  dB)  of  a helical 
antenna  on  an  infinite  ground  plane  radiating  in 
the  axial  mode.  The  observation  plane  is  defined 
by  <)>=90°,  270°  and  0“<e<180°. 


d =0.002X 

a =13.965° 

C = 1.0053  X 
N = 2 


Figure  3-17(d).  The  calculated  E,j  pattern  (not  in  dB)  of  a helical 
antenna  on  a finite  ground  plane  radiating  in  the 
axial  mode.  The  observation  plane  is  defined  as 


d = o.oo2X 

a = 9.0341® 

C = 0.1  ttX 
N = I 

Figure  3-18(a).  The  calculated  Ey  pattern  (not  in  dB)  of  a heliCdl 
antenna  on  an  infinite  ground  plane  radiating  in 
the  normal  mode.  The  observation  plane  is  defined 
by  4>=90°,  270°  and  0°<e<180°. 
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d =0.002X 
a = 9.034r 
C = 0.1  7tX 
N = I 

Figure  3-18(b).  The  calculated  Ee  pattern  (not  in  dB)  of  a helical 
antenna  on  a finite  ground  plane  radiating  in  the 
nonnal  mode.  The  observation  plane  is  defined 
as  in  (a). 


rirjure  3-19(a).  The  calculated  input  impedance  of  a 5-turn 
helix  approximated  by  hexagons  on  an 
infinite  ground  plane  as  a function  of 
frequency. 
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Figure  3-19(b).  The  calculated  input  impedance  of  a b-turn 
helix  approximated  by  hexagons  on  a finite 
ground  plane  as  a function  of  frequency. 


liguro  3-20.  The  Measured  in[)ut  mpodance  of  a b-turn 
helix  on  a ground  plane  as  a function  of 
frequency  is  obtained  by  Kraus. 


IV.  SUMMARY  AND  CONCLUSION 

Antennas  on  satellites  assume  many  different  confiourations.  Many 
of  these  configurations  are  wire  type  antennas  or  can  be  modeled  by 
wires.  In  addition  to  modeling  the  antennas  themselves,  the  satellite 
body  itself  must  also  be  modeled. 

The  work  presented  in  this  final  report  has  shown  how  the  modeling 
of  both  the  antennas  and  satellite  body  can  best  be  accomplished.  The 
ten'is  of  interest  that  can  be  accurately  calculated  include  impedance, 
radiation  patterns,  coupling  between  antennas,  response  of  an  antenna  to 
another  when  the  response  antenna  is  loaded  with  a frequency  sensitive 
device  (filter)  etc. 

The  computer  program  developed  in  connection  with  the  work  discussed 
herein  is  sufficiently  flexible  that  it  can  be  used  in  the  design  as 
well  as  analysis  of  satellite  antenna  systems.  Lor  evample,  the 
performance  of  an  antenna  in  the  presence  of  up  to  four  other  antennas 
of  arbitrary  wire  geometry  (e.o.,  lielix)  can  be  calculated  as  a function 
of  frequency.  The  ultimate  use  made  of  such  software  depends  in  large 
part  upon  the  degree  to  which  the  program  is  made  user  orientated.  One 
possibility  for  the  future  is  to  make  the  program  interactive  through 
a visual  (i.e.,  CRT)  display. 
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